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In an effort to understand the impact of nanostructuring on the magnetocaloric effect, we have 
grown and studied gadohnium in MgO/W(50 A)/[Gd(400 A)/W(50 A)]8 heterostructures. The 
entropy change associated with the second order magnetic phase transition was determined from 
the isothermal magnetization for numerous temperatures and the appropriate Maxwell relation. The 
entropy change peaks at a temperature of 284 K with a value of approximately 3.4 J/kg-K for a 
0-30 kOe field change; the full width at half max of the entropy change peak is about 70 K, which is 
significantly wider than that of bulk Gd under similar conditions. The relative cooling power of this 
nanoscale system is about 240 J/kg, somewhat lower than that of bulk Gd (410 J/kg). An iterative 
Kovel-Fisher method was used to determine the critical exponents governing the phase transition 
to he (3 — 0.51, and 7 = 1.75. Along with a suppressed Curie temperature relative to the bulk, 
the fact that the convergent value of 7 is that predicted by the 2-D Ising model may suggest that 
finite size effects play an important role in this system. Together, these observations suggest that 
nanostructuring may be a promising route to tailoring the magnetocaloric response of materials. 



INTRODUCTION 

Materials exhibiting the giant magnetocaloric effect 
(MCE) have demonstrated potential for advancing mag- 
netic refrigeration, an energy-efficient and environmen- 
tally friendly alternative to conventional refrigeration 
[U [2] . Historically, magnetic refrigerants have been used 
only for very specific applications at very low temper- 
atures, and in most cases, the entropy changes were 
small p . In fact, the largest values of the MCE were 
encountered in the vicinity of a second-order magnetic 
transition in Gd [4j. The discovery of the giant MCE 
in Gd5(Sii_ajGeaj)4 has opened the door for highly effi- 
cient magnetic refrigeration near room temperature [5j|6]. 
However, with the currently available magnetic materi- 
als, this high efficiency is only realized in high magnetic 
fields. Nanostructuring is a promising route to perturb 
properties, and may lead to novel and advantageous mag- 
netocaloric properties using existing materials [7 . As Gd 
is traditionally the standard by which all magnetocaloric 
materials are measured, we have begun investigating the 
impact of nanostructuring using Gd thin films, specifi- 
cally in Gd/W thin film heterostructures. 



EXPERIMENTAL 

Magnetron sputtering was used with ultra high 
purity (UHP, 99.999%) argon gas to deposit 
W(50 A)/[Gd(400 A)/W(50 A)]8 multilayers onto 
MgO (100) substrates at ambient temperature in an 
all-stainless system whose base pressure is 20 nTorr. 
Prior to deposition, the MgO substrate surface was 
modestly milled with 50 W rf biasing in 10 mTorr UHP 
Ar for ten minutes, and the sputtering targets, 99.95% 
W and 99.95% (rare earth equivalent) Gd, were pre- 



sputtered under deposition conditions for ten minutes. 
In addition to cleaning the target surface, presputtering 
the Gd serves to getter background gases from the 
chamber, causing the base pressure to fall by an order 
of magnitude. The deposition parameters were 3 mTorr 
UHP Ar flowing at 20 standard cubic centimeters per 
minute (SCCM) and 100 W of dc power for W, 100 W rf 
for Gd. Rate calibrations of samples grown under these 
conditions were performed on thicker films using x-ray 
reflectivity; the deposition rates were 0.71 A/s for W 
and 1.20 A/s for Gd. The sample was in situ annealed 
at 600 °C for one hour after deposition was complete, 
which is known to improve interface quality [8 . X-ray 
diffraction for 29 in the range 20-70° reveals and strong 
Gd (1010) and (0001) peaks (higher order peaks of these 
are also present). A simple Scherrer grain size analysis 
indicates that the Gd layers are structurally coherent 
throughout their thickness. 

Magnetic measurements were performed using a Quan- 
tum Design Physical Property Measurement System. 
The magnetization isotherms were measured in the range 
of 0-3 T for temperatures of 260-320 K in steps of 10 K. 
The change in the entropy (/S) of a magnetic material in 
an applied magnetic field {H) is related to the change in 
magnetization (M) with respect to the temperature (T) 
through the thermodynamic Maxwell relation: 

fdS{T,H)\ _ fdM(T,H)\ 
\ OH )^-\ dT 

For data taken at discrete field and temperature intervals, 
the change in magnetic entropy, ASm^ due to an applied 
field from to Hq can be approximated as: 

J-i-\-l — -Li 

where fio is the permeability of free space. Figure [l] 
shows the magnetic entropy change measured for the 
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FIG. 1: (Color online) Magnetic entropy change of 
W(50A)/[Gd(400A) /W(50A)]8 . 



W(50 A)/[Gd(400 A)/W(50 A)]8 multilayer (sample area 
was 1 cm^). The ASm peak for low fields is around 
284 K, consistent with Tc as determined by the Kovel- 
Fisher method below [9]. The magnitude of the peak is 
about 3.4 J/kg-K, about one third of the value for bulk 
Gd [10^. The shift of the peak with maximum field is 
approximately linear in the field range studied, with a 
slope of about 0.16 K/kOe. 

Perhaps the most striking feature of the temperature 
dependence of ASm is that its temperature full width at 
half max TpwHM is nearly double that of bulk Gd under 
similar field conditions. In fact, the TpwHM observed 
for AH = 30 kOe exceeds that of bulk Gd for AH = 
50 kOe. While this is compensated for by a reduction 
of the peak value, the relative cooling power cal- 

culated as ASm X TpwHM [in] is 240 J/kg, which is on 
the order of that for bulk Gd (410 J/kg). This may indi- 
cate that nanostructuring is a potential route for develop- 
ing magnetic refrigerants with large useful temperature 
ranges. 

An iterative Kouvel-Fisher method was used to deter- 
mine the critical exponents and Curie temperature of the 
system This approach has been quite successful for 
analyzing amorphous magnetic materials [121 US]- The 
data were initially analyzed on a modified Arrott-Noakes 
plot, MV/5o {H/MY^^% where /^^ = 2/5 and 7o = 4/3 
[M]. The high field portions {H > 1 T) of each isother- 
mal data set were extrapolated to determine the M^-^ 
and (H/M)^-^^ intercepts with second order polynomi- 
als [15], which allows us to determine the spontaneous 
magnetization (Mg) and the inverse of the initial suscep- 
tibility (Xo^)^ respectively. The slopes of the functions 
Y{T) = Ms/{dMs/dT) and X(T) = Xo V(^Xo V^^) 
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FIG. 2: (Color online) Final Y(T) (blue, squares) and X(T) 
(red, circles) used to determine the convergent critical expo- 
nents via the Kovel-Fisher method. The slopes of Y and X 
are respectively and 7~^. The inset shows the evolution 
of the critical exponents through four iterations. 



near Tc are taken as 1//3 and I/7, respectively. The ini- 
tial exponents, Po = 2/5 and 7^ = 4/3, led to nonlinear 
Y{T) and X(T); second order polynomials fit these func- 
tions well, and were used to estimate Pi and 71 near Tc. 
These exponents were then used to create a new modified 
Arrott-Noakes plot (M^/^^ vs (i^/M)^/^i), thus begin- 
ning another iteration of this process. The data were ana- 
lyzed in this fashion for four iterations at which point the 
exponents were unchanging {Ps — and 73 = 74). The 
convergent exponents were determined to be P = 0.51, 
and 7 = 1.75. Figure |2] shows that the final Y{T) and 
X(T) obtained by following this procedure are linear on 
either side of the second order phase transition; the inset 
shows the values of the two exponents after each itera- 
tion. Interestingly, the Curie temperature was relatively 
insensitive to this iterative process, ranging from 282 K 
to 284 K for each iteration. This Curie temperature is 
reduced from that of bulk Gd, which is consistent with 
finite size effects in Gd thin films [16] . Figure [3] shows 
the final modified Arrott-Noakes plot using the exponents 
p = 0.51, and 7 = 1.75. 



DISCUSSION 

The origin of the enhanced TpwHM of the entropy 
change is currently under investigation. The critical ex- 
ponent 7 determined above to be 1.75 is consistent with 
the two-dimensional Ising model, potentially suggesting 
the origin may be related to a change in dimensionality 
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FIG. 3: (Color online) Modified Arrott Plot using the critical 
exponents converged upon after four iterations of the Kouvel- 
Fisher method: /3 = 0.51,7 = 1.75. 



[17]. The suppressed Curie temperature of the Gd does 
indicate conclusively that finite size effects are playing a 
significant role in the system, further lending some cre- 
dence to this possibility. Another potential origin of the 
enhanced TpwHM^ as well as the reduced entropy change 
peak value, is a distribution of Curie temperatures within 
the Gd. Such a situation could arise, for instance, if the 
Curie temperature were suppressed near the Gd/W in- 
terfaces. Indeed, the interfacial interaction between Fe 
and Gd in Fe/Gd superlattices causes the Gd moments 
near the interfaces to exhibit ferromagnetic order above 
the bulk Curie temperature and with a magnetic moment 
exceeding that of bulk Gd [18 . It is therefore conceivable 
that the Curie temperature would be suppressed at the 
Gd-W interface. An investigation of this sort will likely 
require depth profiling by polarized neutron reflectome- 
try [HllO]. 

CONCLUSION 

We have investigated Gd/W multilayers in the context 
of the magnetocaloric effect. It is clear that nanostruc- 
turing Gd significantly impacts the behavior of this ma- 
terial. Relative to bulk Gd, 400 A thick Gd films have 
a reduced entropy change peak value and enhanced en- 
tropy change full width at half maximum. A reduced 
Curie temperature and susceptibility exponent 7 of 1.75 
suggest that finite size effects are impacting the transi- 
tion, and may be related to the departure of the magnetic 



entropy change behavior in thin films from that of bulk 
material. 
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